DESCRIPTION 

METHOD FOR EVALUATING THE DYNAMIC PERSPECTIVE DISTORTION 
OF A TRANSPARENT BODY AND METHOD FOR SUPPORTING THE 
DESIGNING OF A THREE- DIMENSIONALITY CURVED SHAPE OF A 

TRANSPARENT BODY 

TECHNICAL FIELD 

The present invention relates to a method for 
evaluating the dynamic perspective distortion of a 
transparent body and a method for supporting the 
designing of a three-dimensionally curved shape of a 
transparent body. 
BACKGROUND ART 

In manufacturing automobiles or the like, a window 
glass having a complicated three-dimensional shape has 
been used from demands on design and aerodynamics etc. in 
recent years. When an object is seen through such window 
glass with a three-dimensionally curved plane, the object 
is sometimes seen to have a distortion. This phenomenon 
is called a perspective distortion phenomenon, and it is 
known that the phenomenon is caused at a non-parallel 
portion or a curved plane portion in a glass sheet . 
Since the perspective distortion is a factor of hindering 
visibility, in particular, in driving an automobile, the 
permissible maximum value of the perspective distortion 
is ruled in J IS (Japanese Industrial Standards) . 

The perspective distortion is evaluated by 
inspecting the optical performance of an intended glass 



sheet, or examining sensually an intended glass sheet. 
In recent years, there is a technique of evaluating the 
perspective distortion based on a model of three- 
dimensional shape of a glass sheet produced by CAD 
(Computer Aided Design) , instead of the conventional 
optical performance test or sensory examination. 

These conventional evaluation methods evaluate the 
perspective distortion in a case that a stationary object 
is seen through a glass sheet (hereinbelow, referred to 
as the static perspective distortion) . Accordingly, it 
is suitable in a case of evaluating the perspective 
distortion of a glass sheet used for a building. However, 
a problem would arise in a case of evaluating the 
perspective distortion of a glass sheet used for a 
movable body such as an automobile. 

Namely, the evaluation of the static perspective 
distortion is as a result of evaluating the perspective 
distortion at a large number of points in a local portion 
of the glass sheet independently, and it is not as a 
result of considering the continuity of the perspective 
distortion between a plurality of points. Accordingly, 
in a glass sheet wherein the perspective distortion of 
each local portion satisfies a predetermined prescribed 
value, there is a case that a change of the perspective 
distortion between adjacent local portions is large. In 
such case, when a scenery outside the movable body is 
seen through such glass sheet, the perspective distortion 



of an image changes largely with the movement of the 
movable body, and undulation is recognized. 

The important inspection item of a movable body such 
as an automobile is, in addition to evaluating the static 
perspective distortion of an object in a stationary state 
when it is seen through a glass sheet, to evaluate the 
perspective distortion of the object seen through the 
glass sheet in a state of being driven (hereinbelow, 
referred to as the dynamic perspective distortion) . 

Further, in recent years, a technique of designing a 
three-dimensional shape of a glass sheet by CAD has 
widely been used. Accordingly, it is desirable to 
optimize the shape of the glass sheet in consideration of 
an evaluation result of the dynamic perspective 
distortion at the stage of designing. 

The present invention has been achieved in 
consideration of the above-mentioned circumstances, and 
it is a primary object to provide a method for evaluating 
the dynamic perspective distortion of an object in a case 
that a scenery outside the movable body is seen through a 
transparent body such as glass. Further, it is secondary 
object of the present invention to provide a method for 
determining a three-dimensional shape of a transparent 
body based on an evaluation result of the dynamic 
perspective distortion. 
DISCLOSURE OF THE INVENTION 

According to the present invention, there is provided 



a method for evaluating the dynamic perspective 
distortion of a transparent body, which comprises a step 
of producing a model of three-dimensionally curved shape 
of a transparent body having a predetermined refractive 
index; a step of determining an eye point at a side of 
the model of three-dimensionally curved shape and a 
virtual evaluation pattern having a plurality of 
evaluation points at the other side of the model of 
three-dimensionally curved shape; a step of observing, 
from the eye point, the virtual evaluation pattern 
through the transparent body, extracting perspective 
evaluation points as images of the evaluation points, 
obtained by observing through the transparent body, in a 
two-dimensional picture image obtained by the observation, 
and obtaining distance values of adjacent perspective 
evaluation points; a step of determining an optional 
value to be a reference value, among these distance 
values, and a step of evaluating the dynamic perspective 
distortion of the transparent body by obtaining ratios of 
the distance values to the reference value. 

In an aspect of the present invention, there is 
provided the method for evaluating the dynamic 
perspective distortion of a transparent body, wherein the 
dynamic perspective distortion of the transparent body is 
evaluated based on the rate of change of the ratios of 
the distance values to the reference value. 

Further, there is provided the method for evaluating 



the dynamic perspective distortion of a transparent body, 
wherein the minimum value among the distance values is 
selected as the reference value, and the dynamic 
perspective distortion of the transparent body is 
evaluated based on the maximum value among the ratios of 
the distance values with respect to the minimum value. 

Further, there is provided the method for evaluating 
the dynamic perspective distortion of a transparent body, 
wherein the virtual evaluation pattern is an orthogonal 
grid pattern. 

Further, there is provided the method for evaluating 
the dynamic perspective distortion of a transparent body, 
wherein the transparent body is at least one selected 
from a glass sheet and a resinous plate. 

Further, there is provided the method for evaluating 
the dynamic perspective distortion of a transparent body, 
wherein the image seen through the model of three- 
dimensionally curved shape of the transparent body is 
animation-displayed . 

Further, according to the present invention, there 
is provided a method for supporting a model of three- 
dimensionally curved shape of a transparent body, which 
comprises a step of producing a model of three- 
dimensional ly curved shape of a transparent body having a 
predetermined refractive index; a step of determining an 
eye point at a side of the model of three-dimensionally 
curved shape and a virtual evaluation pattern having a 



plurality of evaluation points at the other side of the 
model of three-dimensionally curved shape; a step of 
observing, from the eye point, the virtual evaluation 
pattern through the transparent body, extracting 
perspective evaluation points as images of the evaluation 
points, obtained by observing through the transparent 
body, in a two-dimensional picture image obtained by the 
observation, and obtaining distance values of adjacent 
perspective evaluation points; a step of determining an 
optional value to be a reference value, among these 
distance values; a step of evaluating the dynamic 
perspective distortion of the transparent body by 
obtaining ratios of the distance values to the reference 
value, and a step of correcting the three-dimensionally 
curved shape of the transparent body according to the 
evaluation. 

Further, in an aspect of the present invention, 
there is provided the method for supporting the designing 
of the three-dimensionally curved shape of a transparent 
body, wherein the dynamic perspective distortion of the 
transparent body is evaluated based on the rate of change 
of the ratios of the distance values to the reference 
value . 

Further, there is provided the method for supporting 
the designing of the three-dimensionally curved shape of 
a transparent body, wherein the minimum value among the 
distance values is selected as the reference value, and 



the dynamic perspective distortion of the transparent 
body is evaluated based on the maximum value among the 
ratios of the distance values with respect to the minimum 
value . 

Further, there is provided the method for supporting 
the designing of the three-dimensionally curved shape of 
a transparent body, wherein the virtual evaluation 
pattern is an orthogonal grid pattern. 

Further, there is provided the method for supporting 
the designing of the three-dimensionally curved shape of 
a transparent body, wherein the transparent body is at 
least one selected from a glass sheet and a resinous 
plate . 

Further, there is provided the method for supporting 
the designing of the three-dimensionally curved shape of 
a transparent body, wherein the image seen through the 
model of three-dimensionally curved shape of the 
transparent body is animation-displayed. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 is a flowchart showing the process of the 
method for evaluating the dynamic perspective distortion 
of a transparent body according to an embodiment of the 
present invention . 

Fig. 2 is a constitution diagram showing an 
evaluation system for carrying out the evaluation method 
shown in Fig. 1. 

Fig. 3 is a constitution diagram showing an 
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evaluation model. 

Fig. 4 is a flowchart showing the process of 
calculation for obtaining a deviation of the position of 
a grid point seen through the transparent body. 

Fig. 5 is a diagram showing the process of 
calculation for obtaining a deviation of the position a 
grid point seen through the transparent body. 

Fig. 6 is a diagram showing an evaluation pattern 
formed freshly according to the evaluation model . 

Fig. 7 is a graph showing grid distance ratios in 
lateral directions . 

Fig. 8A is a graph showing a result of calculation 
of dynamic distortion values. 

Fig. 8B is a graph showing a result of measurement 
of the dynamic perspective distortion in a sensory 
evaluation test. 

Fig. 9 is a graph showing the relation between the 
dynamic distortion value of the transparent body and the 
intensity of the dynamic distortion. 
BEST MODE FOR CARRYING OUT THE INVENTION 

In the following, some embodiments of the present 
invention will be described with reference to the drawing 

Fig. 1 is a flowchart showing the method of 
evaluation of the dynamic perspective distortion 
according to an embodiment of the present invention. Fig 
2 is a block diagram showing an evaluation system for 
carrying out the method for evaluating the dynamic 



perspective distortion shown in Fig. 1. Fig. 3 shows an 
evaluation model. 

First, a model of three-dimensionally curved shape 
of a transparent body (such as a glass sheet or a 
resinous plate) is produced by reading out the data of 
shape by means of an input device, a memory device or the 
like (Step SI) . The method for inputting the data of 
shape is well known, and CAD data 5 may be inputted into 
a computer 9 with a magnetic tape or the like by means of 
a magnetic reproducing device 6, or data obtained by 
plotting a designed drawing 7 with a digitizer 8 may be 
inputted into the computer 9, as shown in Fig. 2. 
Further, a known technique can be used for producing the 
three-dimensionally curved shape of the transparent body. 
For example, three-dimensional coordinate data of a large 
number of points on a plan view and a side view, called a 
Mylar diagram, of the transparent body are inputted into 
the computer 9 . 

Then, a central processing unit 17, which has 
received input data, produces grid-like three-dimensional 
spline curves passing a train of points; divides a curved 
plane into rectangular patches defined by the spline 
curves as boundaries, and produces a dual third 
parametric curved plane like a Coons plane which is 
expressed by parameters along each side of patches. And, 
it renders the curved plane to be the shape model of an 
inner surface of a plate-like transparent body, and 
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produces the shape model of an outer surface by 
offsetting the shape model of the inner surface by the 
thickness of the transparent body. The finally 
determined shape model is obtainable by combining the 
shape models of the inner and outer surfaces. 

The data of the obtained shape model are stored in a 
shape data file 18. The above-mentioned series of works 
for producing the model of three-dimensionally curved 
shape are carried out by the central processing unit 17 
which reads out for execution a program stored as a shape 
model producing means 11 in a program section 10 of the 
computer 9 . 

The central processing unit 17 determines an eye 
point EP (Step S2) and produces virtual grid board (Step 
S3) based on data in a condition data file 19 in which 
determination conditions previously inputted by a key- 
input means 4 are stored. In this embodiment, a virtual 
grid board 2 (see Fig. 3) is assumed for a virtual 
evaluation pattern, and grid points (PS1, PS2 and PS3 ) 
are assumed as evaluation points. However, another 
evaluation pattern having regularly arranged points in 
lateral and vertical directions may be used. Evaluation 
points seen from the eye point EP become grid points (PRl, 
PR2 and PR3 (hereinbelow, referred to as perspective 
evaluation point)) due to the refraction index of the 
transparent body. 

Further, based similarly on the data in the 
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condition data file 19, the fitting angle of the shape 
model of the transparent body; the distances between the 
eye point EP and the shape model 1 of the transparent 
body and between the eye point and the virtual grid board 
2, and the grid distance of the virtual grid board 2 are 
determined. As described above, all conditions before 
the initiation of the measurement are determined, and 
obtained arrangement data are stored in an arrangement 
data file 20 (Step S4) . 

A series of these works for determining evaluation 
conditions are carried out by the central processing unit 
17 which reads out for execution a program stored as an 
evaluation condition determining means 12 in the program 
section 10 of the computer 9. 

Then, the computer 9 calls an optical path 
tracing/calculating means 13 in the program section 10 to 
define a virtual line with respect to the shape model 1 
of the transparent body by using the arrangement data on 
the shape model 1 of the transparent body, the eye point 
EP and the virtual grid board 2 (Step S5) , and to conduct 
calculation for obtaining a perspective distortion due to 
the refraction (a deviation of the position of a grid 
point observed from the eye point EP) (Step S6) ♦ The 
calculation of the deviation of the position of a grid 
point is, in principle, conducted by utilizing the fact 
that the propagating direction of a light beam is 
different due to the refraction effect depending on the 
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presence or absence of the shape model 1 of the 
transparent body. 

Fig. 4 shows concretely the flow of the calculation. 
Further, the tracing of the optical path is explained in 
detail by using Fig. 5. First, an optional grid point PO 
is determined on the virtual grid board 2 (Step SSI) . 
Then, a vector VRO in the propagating direction of a 
virtual light beam 3 directing from the eye point EP 
toward PO is obtained (Step SS2). Then, the intersection 
point PI between a linear line passing through the eye 
point EP and having the same direction as the vector VRO 
and an inner surface SU1 of the shape model 1 of the 
transparent body is obtained (Step SS3), and a vector VE1 
along the normal line of the inner surface of the shape 
model 1 of the transparent body at the intersection point 
Pi is obtained (Step SS4) . 

Further, based on the vector VRO and the vector VEl 
and according to the law of refraction, a vector VRl 
which indicates the propagating path after the virtual 
light beam 3 is refracted at the outer surface SU2 of the 
shape model 1 of the transparent body, is obtained (Step 
SS5) . Namely, when the incident angle (with respect to 
the normal line) of the virtual light beam 3 from the eye 
point EP to the inner surface SUl, is assumed i, the 
refractive angle (with respect to the normal line) of the 
virtual light beam 3 at the inner surface SUl is r and 
the refractive index of the transparent body to air is n, 



sin i/sin r=n is established, and accordingly, the 
refractive angle r can be obtained from the known vector 
VRO, vector VE1 and the refractive index n. Thus, the 
vector VR1 is obtainable based on these. 

Then, the intersection point P2 between a linear 
line passing through the intersection point PI and having 
the same direction as the vector VR1 and an outer surface 
SU2 of the shape model 1 of the transparent body is 
obtained (Step SS6) , and a vector VE2 along the normal 
line on the outer surface SU2 of the shape model 1 of the 
transparent body at the intersection point P2 is obtained 
(Step SS7). Further, based on the vector VR1 and the 
vector VE2 and according to the above-mentioned law of 
refraction, a vector vR2 which indicates the propagating 
path after the virtual light beam 3 is refracted at the 
outer surface SU2 of the shape model 1 of the transparent 
body, is obtained (Step SS8) . Then, the intersection 
point P3 between a linear line passing through the 
intersection point P2 and having the same direction as 
the vector VR2 and the virtual grid board 2, is obtained 
(Step SS9) . 

After this, the position of Pi is shifted so that P3 
approaches PO, and P3 is again obtained according to the 
above-mentioned Steps SSI to SS9 . Namely, Pi is moved by 
a distance shorter than a line segment P3P0 in the 
direction of a vector P3P0 to determine a fresh Pi 
(hereinbelow, referred to as Pli) , so that P3 is again 
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obtained. As a result, when the line segment P3P0 is 
shorter than a previously determined length (hereinbelow, 
referred to as P3 = P0), Pli is replaced by PI, and VR2 is 
obtained. If P3^P0 is not established, a fresh Pi 
(=P1 2 , Pin (n: an arbitrary natural number)) is 

determined to repeat the above-mentioned steps until P3^ 
PO can be established (Step SS10). 

As described above, the position of the intersection 
point (virtual point) PI can be in approximation the 
position where the virtual point PO can actually be 
observed from the eye point EP. 

With respect to all the other grid points, 
intersection points corresponding to respective grid 
points are obtainable by executing Steps SSI to SS10 
(Step S7) . By connecting these intersection points, a 
fresh shape of the grids can be constructed (Step S8) . 
The data of the intersection points are stored in a 
result data file 21. 

Then, the dynamic distortion value as a parameter 
for evaluating the dynamic perspective distortion and the 
intensity of the dynamic distortion value are obtained. 
In this embodiment, the program for obtaining each 
parameter is stored in a result data evaluation means 14 
in the program section 10. First, based on the above- 
mentioned result data, the narrowest grid distance in a 
vertical direction among grid distances (between the 
intersection points) which are seen from the transparent 



body, is calculated. For example, explanation is made 
assuming that perspective evaluation points PRl, PR2 and 
PR3 as intersection points corresponding to three grid 
points PS1, PS2 and PS3 in Fig. 3, have been obtained. 
In this case, the distance between the perspective 
evaluation points PRl and PR2 is the grid distance in a 
vertical direction. The distance can be obtained from 
the coordinate values of two points. By obtaining all 
grid distances in the vertical direction as described, 
the minimum value among them can be found (Step S9) . 

Fig. 6 shows an evaluation pattern formed freshly 
from the evaluation model wherein dl to dn indicate grid 
distances in a vertical direction. In Fig. 6, the 
minimum value of the grid distance is indicated by d4 . 
By dividing each of the other grid distances by the 
minimum value, grid distance ratios di/d4 (i=l to n, n: 
the number of grids) can be obtained (Step S10) . The 
maximum value of the grid distance ratios is determined 
to be a dynamic distortion value (Step Sll) . 

Then, a dynamic distortion value is also obtained in 
a lateral direction in the same manner. First, a grid 
distance in a lateral direction is obtained. In the 
example of Fig. 3, the distance between perspective 
evaluation points PRl and PR3 is the grid distance in the 
lateral direction. By obtaining all grid distances in 
the lateral direction as described, the minimum value 
among them can be found. By dividing each of the other 
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grid distances by the minimum value, grid distance ratios 
can be obtained. The maximum value among them is 
determined to be a dynamic distortion value (Step S12) . 

In this embodiment, the minimum value of grid 
distances is determined to be the reference value and 
each grid distance ratio can be obtained by dividing each 
of the other grid distances in vertical and lateral 
directions by the reference value. The maximum value 
among them is defined as the dynamic distortion value. 
However, the dynamic distortion value is a relative value 
based on the ratio of the reference value which is an 
optional one among the grid distances, to the other grid 
distance. Accordingly, the dynamic distortion value may 
be determined irrespective of this embodiment. As in 
this embodiment, when the reference value is determined 
as the minimum value of the grid distances and the 
maximum value is determined in the relation of the other 
grid distances to the minimum value, the calculation in 
the determination of the minimum value can quickly and 
easily be performed. Further, the difference between an 
expansively visible portion and a diminutively visible 
portion due to the perspective distortion can easily be 
expressed numerically. 

Then, the rate of change of the grid distance ratios 
in both the vertical and lateral directions (two- 
dimensional direction) is calculated (Step S13). For 
example, the grid distance ratios are placed in the order 
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from left to right in the direction of a linear line Ll 
which traverses laterally the virtual grid board 2 in Fig. 
3. Fig. 7 is a graph showing such state which is 
obtainable by selecting a region wherein the gradient 
seems to be the largest, and applying a least squares 
method or a suitable method to the grid distance ratio of 
this region. The rate of change of the grid distance 
ratios is calculated as the gradient based on the graph, 
and the maximum value of the gradient is determined as 
the intensity of the dynamic distortion (Step S14) . As 
the rate of change is larger, the difference of the 
perspective distortion between that portion and 
surroundings is larger. Even though there are some 
dynamic distortion values (the maximal grid distance 
ratio values) , flickering in view of an object seen 
through the transparent body from a movable body is apt 
to be realized if the rate of change is larger. By the 
introduction of this rate of change, the continuity of 
the perspective distortion is taken into consideration 
whereby the evaluation of the dynamic perspective 
distortion becomes possible. 

Thus, the processes for obtaining the dynamic 
distortion value and the intensity of the dynamic 
distortion value are finished about a single shape model. 
The obtainable result is preserved in the result data 
file 21 by means of the result data evaluation means 14. 

It is desirable that the data of grid distance 
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ratios, the dynamic distortion values, the rates of 
change of the grid distance ratios and so on, which are 
obtained respectively, are displayed by a result display 
means 22 or printed by a printer 23 occasionally in the 
form of a distribution map of the grid distance ratios or 
a graph of the rates of change of the grid distance 
ratios, although the flowchart in Fig. 1 does not show 
specifically. 

Then, the computer calls a model shape determining 
means 15 in the program section 10. In this program, the 
obtained dynamic distortion value is compared with a 
previously determined threshold value (Step S15) . When 
it is lower than the threshold value, a judgement of 
acceptance is made whereby the operation of judgement is 
finished. When it is higher than the threshold value, a 
judgement of rejection is made whereby the operation of 
shape adjustment is carried out on CAD (Step SI 6) . In 
the adjustment, the portion of the shape model which 
corresponds to a portion where the grid distance becomes 
maximal, is displayed, and it is requested to change the 
surface curvature and the trim shape of that portion. 
After the adjustment, the above-mentioned data processing 
is carried out to obtain the dynamic distortion value 
(Step S17) . Simulation and adjustment are repeated until 
a value lower than the threshold value can be obtained. 
Thus, the optimum shape can be obtained. 

Finally, the computer 9 calls an animation producing 
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means 16 in the program section 10. A scenery seen 
through the accepted shape model of the transparent body 
from, for example, a cruising automobile is displayed in 
the form of an animation image in a graphic display 2 4 in 
order to confirm visually that the adjustment has been 
completed. The animation display is conducted by 
utilizing the result data stored in the result data file 
21. The operation for the shape model before the 
adjustment is also conducted similarly. Namely, a 
scenery seen through the transparent body from, for 
example, a cruising automobile is displayed in the form 
of an animation image such as CG (Computer Graphics) or 
the like, in the graphic display 24. As a result, a user 
can confirm effectively adjustment effect by comparing 
the animation image after the adjustment with that before 
the ad j us tmen t . 

In the following, explanation will be made as to a 
method for determining the threshold value in the 
judgement of shape which is carried out in this 
embodiment. Shape models were prepared using a large 
number of actual transparent bodies, and the dynamic 
distortion values were calculated on the shape models 
respectively (Fig. 8A) . At the same time, sensory 
evaluation tests for 5-step evaluation were conducted to 
a plurality of tested persons to examine the feeling of 
flickering in a case that a moving object was actually 
seen through these transparent bodies (Fig. 8B) . 
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In both graphs, the abscissa indicates samples A to 
I and the ordinate indicates the dynamic distortion value. 
As clear from the comparison of these graphs, it could be 
confirmed that the dynamic distortion value has extremely 
high correlation to the evaluation by actual observation 
by eye, namely, the value was an index which sufficiently 
reflected the dynamic perspective distortion. Then, 
based on the fact that the sample H and the sample I were 
judged as no good in the sensory evaluation tests, the 
threshold value of the dynamic distortion value was 
determined to be 1.6. 

In the measurement of the intensity of the dynamic 
distortion (the maximum gradient of the grid distance 
ratios) , in addition to the dynamic distortion value (the 
maximum value of the grid distance ratios) , on the large 
number of transparent bodies, there was found 
distribution in a region CD to a region ® in Fig. 9. 
From the visual confirmation by the above-mentioned 
animation and the actual evaluation, the range of a 
region A in which flickering was hardly realized was 
determined to be the reference of judgement. By 
conducting the double judgement in terms of the dynamic 
distortion value and the intensity of the dynamic 
distortion, it is possible to obtain a highly accurate 
evaluation in a step of designing. 
INDUSTRIAL APPLICABILITY 

According to the present invention, it is possible 
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to quantify the dynamic perspective distortion (the 
dynamic distortion value and the intensity of the dynamic 
distortion) of a transparent body, which has not been 
used for the evaluation, and evaluation can be conducted 
5 with higher reliability than conventional technique to 
the dynamic perspective distortion which is caused when 
outward things are seen through a transparent body such 
as glass, resin or the like from a movable body such as 
an automobile. Further, the feedback to the 

10 manufacturing process can be easily with the improvement 
of accuracy in designing; the yield percentage and the 
quality can be improved; cost for unnecessary manufacture 
of mold or modification of mold can be reduced to reduce 
cost, and the flexibility of designing can be expanded. 

15 Accordingly, the present invention is in particular 
effective for large item small scale production. 

The entire disclosure of Japanese Patent Application 
No. 2000-193035 filed on June 27, 2000 including 
specification, claims, drawings and summary are 

20 incorporated herein by reference in its entirety 
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direction 
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Fig. 4 (attached sheet) 



551 Determining PO 

552 Obtaining VRO directing from EP toward PO 

553 Obtaining an intersection point Pi at the intersection 
of VRO to an inner surface SU1 of the shape model 

554 Obtaining a normal line vector VEl on the inner 
surface SU1 at PI 

555 Obtaining a refraction vector VRl from VRO and VEl 

556 Obtaining an intersection point P2 at the intersection 
of VRl to an outer surface SU2 of the shape model 

557 Obtaining a normal line vector VE2 on the outer 
surface SU2 at P2 

558 Obtaining a refraction vector VR2 from VRl and VE2 

559 Obtaining an intersection point P3 at the 
intersection of the line passing through VR2 to a 
grid broad 2 

SS10 Moving repeatedly Pi to find VR2 so as to approach 
P3^=P0 r and repeating the above operation to find 
PI so as to approach P3^P0 



